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Abstract 



In this report we discuss the impact of polarized foregrounds on a future CMBPol satellite mission. We 
review our current knowledge of Galactic polarized emission at microwave frequencies, including synchrotron 
and thermal dust emission. We use existing data and our understanding of the physical behavior of the sources 
of foreground emission to generate sky templates, and start to assess how well primordial gravitational wave 
signals can be separated from foreground contaminants for a CMBPol mission. At the estimated foreground 
minimum of ~ 100 GHz, the polarized foregrounds are expected to be lower than a primordial polarization 
signal with tensor-to-scalar ratio r = 0.01, in a small patch (~ 1%) of the sky known to have low Galactic 
emission. Over 75% of the sky we expect the foreground amplitude to exceed the primordial signal by about 
a factor of eight at the foreground minimum and on scales of two degrees. Only on the largest scales does the 
polarized foreground amplitude exceed the primordial signal by a larger factor of about 20. The prospects 
for detecting an r — 0.01 signal including degree-scale measurements appear promising, with 5ov ~ 0.003 
forecast from multiple methods. A mission that observes a range of scales offers better prospects from 
the foregrounds perspective than one targeting only the lowest few multipoles. We begin to explore how 
optimizing the composition of frequency channels in the focal plane can maximize our ability to perform 
component separation, with a range of typically 40 ^ v ;$ 300 GHz preferred for ten channels. Foreground 
cleaning methods are already in place to tackle a CMBPol mission data set, and further investigation of the 
optimization and detectability of the primordial signal will be useful for mission design. 
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1 Introduction 



Measurements of the Cosmic Microwave Background (CMB) temperat ure anisotropy have led to the es 



tablishment of a standard ACDM cosmological model for the univers e (Miller et al 



Netterfield et al 



2002; 



Halverson et al. 2002: Pearson et al. 2003: Scott et al 



200. 



3: 



1999; 



Lee et al 



Spergel et al 



200. 
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2001 



2007 



Dunklev et al.l l2008al) . The pola rization anisotropy, two orders of m agnitude smaller, was measured for 



the first time in 2002 by DASI |Kovac et all l200a iLeitch et all l2002h and confirmed by t h e ground and 



balloon-based experiments CBI, CAPMAP, Boomerang, and QuAD (IReadhead et alj l2004t iBarkats et al 



2005; Montrov et al 



2006. P. A de et alJl2 007: Prv ke et alJl2008h . The first all-sky observations by the WMAP 



satellite (jPage et alJl2007t iHinshaw et al.ll2008h . with a measurement of the power in the curl- free polarized 
'E-mode', have provided a cross-check of the cosmological model, led to improved constraint s on cosmologica l 
parameters, and a m easurement of the optical depth to the reionization of the universe (jPage et al.l 12007 : 



Hinshaw et al 



2008h . 



The next challenge of CMB observation is to test the inflationary scenario for the early universe, by looking 
for the signature of primordial gravitational waves. They are predicted to leave a distinct divergence-free, 
or 'B-mode', pa ttern in the large-scale CMB polarization anisotropy that may be observable with a future 
experiment (see lBaumann et al. f (|2008bh for details). This signal is already constrained to be over an order of 
magnitude smaller t han the observed E-mode signal arising from scalar fluctu ations, from ground and space- 
based observations ( Page et al. 2007 : Hinshaw et al. 20081 : Prvke et al. 2008 ). The Planck satellite, due for 
launch 2 009, is projected to reach limits on the ratio of power in primordial tensor to scalar fluctuations r, of 
r ~ 0.1 (jThe Planck Collaborationll2006l ) . while ground and balloon-based experiments currently observing or 
being co mmissioned a r e desi g ned to reach below r = 0.1 . Thes e include EBEX, C lover, QUIET, and SPIDER 
(see e.g., bxlev et al.l |2004h ; iTavlorj (|2006h ; ISamtlebenl |2008h ; ICrill et all |2008h for descriptions). 

This document is part of a study focusing on the capabilities of a future satellite mission targeting limits 
of r = 0.01 or lower, known as CMBPol. Here we address the issue of polarized foreground emission from the 
Galaxy, which on average dominates the primordial B-mode over the whole sky. Establishing how well we can 
realistically hope to extract the primord ial signal, and how t h is influences mission design, has b een a d dress ed 
in a number of previous studies (e.g., lAmarie et aD |2005h ; IVerde et al.l |2006h ; bock et all h00& l2008h). 



In com panion documents the theoretical c ase for inflat i on an d its predictions are presented (Bauman n et al 



2008b), with weak lensing effects studied in lSmith et al.l (12008!). and r eionizat ion prospects in Zaldar riaga et al 
(j2008h . The Galactic science case is presented in Fraisse et al. ( 2008h . and Baumann et al. |2008a ) provide a 
summary. 

The document is structured as follows. In fj2]we review our current understanding of polarized foreground 
emission, and discuss predictions and models for the emission, with estimates of large-scale polarization maps 
and power spectra. In Sj3] we describe various methods that are used to estimate polarized CMB maps, and 
in fJH apply these and Fisher matrix methods to begin to forecast limits for a specific template mission, and 
test optimal allocation of channels among frequencies. We also describe future tests, and summarize in S}5] 



2 Knowledge of polarized foregrounds 

In this section we focus on observations of polarized foregrounds relevant to the microwave regime, concen- 
trating on large-scale diffuse emission. The discussion is organized by physical emission mechanism, and 
outlines how our understanding of the physics of the interstellar medium is used to construct two-dimensional 
emission templates on the sky. This is an inherently simpler problem than a three-dimensional model. Issues 
that remain open for investigation include: possible spatial and frequency variation of spectral indices, the 
number of parameters required for a sufficiently sophisticated model to describe and subtract foregrounds at 
the level required if r < 0.01, and methods for assessing the validity of the chosen model. 
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Hinshaw et al. 2008 



Figure 1: The WMAP measurements of polarization at 23 GHz map the large-scale distribution of 
polarized synchrotron emission (Mollweide projection in Galactic coordinates). White lines indicate 
the polarization angle 7^ and amplitude, and are oriented perpendicular to the magnetic field. 

2.1 Synchrotron emission 

Synchrotron emission results from the acceleration of cosmic-ray electrons in the magnetic field of the Galaxy. 
For a power-law distribution of electron energies, N(E) oc E~ p , propagating in a uniform magnetic field, the 
resulting emission is partially polarized with fractional linear polarization 



aligned perpendicular to the magnetic field (jRvbicki fc Lightmanl 119791 ). The frequency dependence of syn- 
chrotron emission is also related to the electron energy distribution, T(v) oc v" , with spectral index 

P = ~ P -^ (2) 

where T is in units of antenna temperature. For spectral index (3 ~ —3 observed at microwave frequencies, 
synchrotron emission could have fractional polarization as high as f s ~ 0.75, although in practice this is 
almost never observed. Line-of-sight and beam averaging effects tend to reduce the observed polarization by 
averaging over regions with different electron energy distribution or magnetic field orientation. At frequencies 
below a few GHz, Faraday rotation will induce additional depolarization. 
The polarized synchrotron emission may be parameterized as 

/ J,(«)+C(S)los(»W 
Q(n) = Tf(n)co S (2 7 (n))(-J 

U{h) = rf(n)sin(2 7 (n))(-J (3) 
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Figure 2: Power spectra for the WMAP 3-year polarization maps are dominated by foreground 
emission. Solid lines show EE spectra and dashed lines show BB spectra outside the Galactic plane; 
the color indicates frequency band. The rise in power for £ > 100 is an artifact of the instrument noise. 
The red dot-dashed line shows the estimated BB foreground power at 60 GHz, using a parameterized 
model fit to the multi-frequency data. Diamonds (EE) and boxes (BB) i ndicate points wit h negative 
values. The cosmological model shown has r = 0.3 and r = 0.09. From lPage et al.l (|2007l ). 



where Q(fi) and U(h) are the Stokes parameters in direction h, Tf(n) is the (polarized) amplitude T p = 
\J Q 2 + U 2 , j(n) is the polarization angle, f3 s (n) is the spectral index, and C(n) parameterizes any spectral 
curvature (flattening or steepening) relative to a pure power law. 

The Wilkinson Microwave Anisotropy Probe (WMAP) provides our best estimate of the synchrotron 
morphology on angular scales of a few degrees o r larger. Figure [1] shows the full-sky WMAP observations 
at 23 GHz (jPage et al.ll2007t iHinshaw et alj |2008). At 4° pixelization, the signal to noise ratio for polarized 
synchrotron emission is better than three over 90% of the sky. The polarization angle has coherent structure 
over large swaths of the sky, creating significant emission at low multipole moment s I. Figure [2] shows the 
power spectra of the WMAP polarization data from 23 to 94 GHz (jPage et al.ll2007[ ). The power spectra are 
dominated by foreground emission with roughly equal power in E- and B-modes, which, averaged over the 
high-latitude sky, are brighter than the CMB polarization even at the minimum near 70-100 GHz. A simple 
parameterization of the foreground emission in lPage et al. ( 2007 ) gives a power spectrum with £(£+l)Ci/2ir cx 
£~°' e , shown in red in FigureO This is consisten t with the power spectrum for synchrotron emission observed 



from radio maps (see e.g. lLa Porta et al 



(200 



We can estimate the fractional polarization of synchrotron emission by comparing the WMAP 23 GHz 
map to an estimate of the unpolarized synchrotron emission. The synchrotron intensity suffers from confusion 
between synchrotron, free-free, and other microwave compone nts, so extracting this signal from observed sky 
maps depends on modeling assumptions. iKoeut et alJ (|2007|) use the WM AP three-year maximum-entropy 
map of synchrotron emission as a tracer of the intensity (jHinshaw et al 1l200i . finding the Galactic plane region 
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Figure 3: Estimated polarization fraction / = P/I at 5° resolution. The polarization P 



\J Q 2 + U 2 is well estimated from the WMAP K band data ( Page et al.1 120071 ) and the intensity 
/, which is less certain, is m odeled excluding (top), and inc luding (bottom), an unpolarized 'anoma- 
lous' dust component (from iMiville-Deschenes et al.l (|2008l )). 



(|6| < 5°) to be largely depolarized, with mean fractional polarization f s = P/I = 0.05. The North Galactic 
spur region at mid-latitudes has fractional polarization /■ 0.3 while the high latitude sky outside the P06 
mask has a very broad distribution with mean f s = 0.15 (jKogut et al.ll2007l K The synchrotron polarization is 
consistent with the magnetic field having rough equ i partition between a large - scale smooth component and a 
smaller scale turbulent component ( Page et al.l2007 ) . Miville-Deschenes et al. ( 20081 ) estimate the polarization 
fraction at 23 GHz for a set of synchrotron intensity maps. One map is derived from WMAP data assuming 
that no anomalous dust emissio n contributes at 23 GHz, and gives a polarization fraction, shown in Figure 
02 similar to the one studied in iKogut et al.l ( 2007 ). An alternative intensity map is derived from WMAP 
data allowing for an anomalous unpolarized emission component to contribute to the total observed intensity. 
This results in a modified estimate for the synchrotron intensity, and an increase in polarization fraction over 
much of the sky, also shown in Figure[3l In this second map, there is a larger polarization fraction in the anti- 
center region compared to the inner galaxy. These estimates for the polarization fraction would be modified if 
the anomalous component (s) had a non- negligible polarization fraction. We note however that despite these 
uncertainties in the polarization fraction, our estimate of the large-scale polarization of synchrotron emission, 
P = yfQ 2 + U 2 , is robust at 23 GHz. 

The WMAP data provide a reasonable signal to noise ratio for the polarized synchrotron morphology on 
angular scales of a few degrees or larger. The spectral index f3 s (h) is less well constrained. Estimates of the 
spectral index for unpolarized synchrotron emission suffer from confusion with other emission components, 
while estimates based on polarization data (where confusion is less severe) suffer from the weaker synchrotron 
signal at frequencies above 23 GHz. Figure 2] shows the observed polarization variance (P 2 ) of the WMAP 
data as a function of frequency for 4° pixelization. The spectra are consistent with a superposition of two 
power-law components, with synchrot ron S » — 3 and d ust (3d ~ 2. The CMB contributes lower power 
relative to the foreground components. IKogut et al. ( 2007 ) find t hat the synchrotron spectral index steepens 



off the plane by a modest amount A/3 S ~ 0.2. iGold et al.l ( 2008f) map the synchrotron spectral index using 
a combined analysis of the temperature and polarization data, finding a trend towards a steepening of the 
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Figure 4: Spectrum of the variance in the WMAP polarization data, computed at each frequency 
from the cross-spectra of one-year maps, to avoid a bias from instrument noise. The spectra are 
consistent with power-law synchrotron and dust emi ssion, with synchro tron (3 S ~ —3. The P06 cut 
excludes 25% of the sky in the Galactic plane. From lKogut et al.1 (120071 ). 



spectral index off the Galactic plane. A single index of S = —3 .3 currently works sufficiently well outside the 
Galactic plane for cleaning the WMAP maps (jPage et al.ll2007J ). 

Similar problems affect estimates of the synchrotron spectral curvature. The spectral index depends on 
the energy spectrum of cosmic ray electrons (Eqn.[2]). Higher energy electrons lose energy more rapidly, which 
can lead to a steepening or break in the synchrotron spectrum. The frequency of the break is not well known, 
but is thought to vary with position on the sky and be in the range 10-100 GHz, corresponding to a break in 
the electron spectrum at the energy where the characteristic cooling time and the escape time are the same. 
To a large extent, efforts to characterize any spectral steepening are limited by confusion with thermal dust 
emission. The effects of steepening for the synchrotron component are largest near the foreground minimum 
near 70-100 GHz, where the amplitudes of synchrotron and thermal dust emission are comparable. New 
measurements of polarized dust emission at frequencies above 100 GHz, from the Planck satellite and other 
experiments, will help break this degeneracy and provide tighter limits on both the dust and synchrotron 
spectra. Low frequency observations from experiments at 5-15 GHz, including the C-Band All Sky Survey 
(C-BASS), will also help determine the frequency dependence of the polarized synchrotron emission. 

Little information is available for synchrotron polarization on angular scales below 1°. Existing mea- 
surements on small scales are heavily affected by Faraday rotation at low frequencies. Figure [5] shows a 
olarization map with angular resolution 36' from the Dominion Radio Astrophysical Observatory at 1.4 GHz 
20061 ). Although some features such as the North Galactic spur at high latitude and the Fan 



Wolleben et al 



region near longitude 150° are recognizable from the polarization map at 23 GHz (Fig[T]), the 1.4 GHz map 
is dominated by Faraday depolarization, whereas the 23 GHz emission has almost no Faraday rotation. Sim- 
ilar depolarization affects interferometric maps at finer angular resolution ( Uvaniker et al.l 2003t Tavlor et al 
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Figure 5: Polarization intensity measured at 1.4 GHz. The prominent Faraday depolarization 
limits use of low- frequ ency maps for foreground modeling, despite the finer angular resolution. From 
Wolleben et al.l (|2006l ). The gray-scale runs from to 450 mK in brightness temperature. 



20031 ). limiting their utility for foreground modeling at millimeter wavelengths. 



An additional component of synchrotron emission that may be constrained with future polarization data 
is anomal ously hard s ynchr o tron emission towards the Ga la ctic center. This is k nown as the WMAP 'Haze', 
found bv iFinkbeineil l|2004l ); bobler & Finkbeinei] (|2008ah ; iBottino et all |2008) in analyses of the WMAP 
data. This does not have a conclusively confirmed origin, but as it is confined to t he Galactic center, it is 
discussed in more detail in the companion CMBPol document on foreground science ( Fraisse et al. 20081 ). 



2.1.1 Modeling considerations 



Given the spatial variation in the observed synchrotron polarization fraction, it is unrealistic to model the 
polarized signal as a synchrotron intensity scaled by a single global polarization fraction, as illustrated in 
Figure [3] This implies that for a given frequency, the emission in every pixel should be modeled with at least 
two parameters for the synchrotron Q and U (or P and 7). In terms of the spectral behavior, th e power law 



approximation seemed to be sufficient for cosmological analyses well outside the plane for WMAP (jPage et al 



20071 : iGold et al.l 120081 ). However, a power law fit in smaller pixels, and a curvature fit either globally or in 
pixels, are likely necessary when pushing the foregrounds down another factor of 10-100. Values for the 
synchrotron index in the range —3.5 f3 a — 2.5 are physically reasonable, with index curvature expected 
to lie in the range —0.5 i£ C 0.5, consistent with current observations. Additional parameters beyond the 
curvature may be required to quantify the frequency dependence, although they are not necessary to fit 
current data. Diffusion of electrons limits the possible spatial variation of the spectral index at small angular 
scales (see e.g.. IStrong et al.l (|2007l )). so models where the index is fit in larger regions may be considered. In 
principle the spectral indices of the observed synchrotron intensity and polarization should be modeled with 
independent parameters, due to possible depolarization effects from the Galactic magnetic field, although a 
difference between the two has not yet been conclusively detected. This is due to both uncertainty in the 
spectral index of the synchrotron intensity, and the instrumental noise levels of the polarization measurements. 



2.2 Dust emission 

Galactic emission in the 100 - 6000 GHz frequency range (A = 50 - 3000 fim) is dominated by thermal 
emission from warm (~ 10 - 100 K) interstellar dust grains. Here 'thermal' refers to emission from ther- 
mal fluctuations in the electric dipole moments of grains. The spectral shape of this emission is generally 
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Figure 6: Thermal dust intensity observed by IRAS and extrapolated to 94 GHz by lFinkbeiner et al- 
ii 19991 ) . This map is commonly used to identify regions of expected low dust emission, with the 
minimum at Galactic coordinates (I, b) ~ (240,-70). Although the dust polarization fraction is 
expected to vary spatially, this map is also useful for identifying regions of potential low dust 
polarization. 



modeled with one or more thermal components modified by a frequency-dependent emission. That is, for a 
single temperature, the intensity is proportional to I(y) ~ v 13 B V (T) 1 where B U (T) is the Planck function at 
frequency v and temperature T . Multiple temperatures and spectral indices (/3) are often needed to model 
the intensity spectrum a t any single point on the sky. Using data from the all-sky surveys of IRAS at 100 
|Neugebauer et al.lll984h . and COBE/DIRBE at 100 and 240 fim (jFixsen et al.lll998h . lFmkbeiner et aD(|l999h 
have modeled this Galactic dust emission using two temperature components, with (T\^) — 9.5 and 16 K; 
/3i,2 = 1-7 and 2.7. Both components are in equilibrium with the interstellar radiation field. They then used 
this model to predict thermal dust emission at microwave frequencies. The resulting model is calibrated from 
COBE/DIRBE and has the angular resolution of IRAS 6'). The predict ed signal, shown in Figure [6l has 
been found consistent with WMAP total intensity data ([Bennett et al.ll2003f ). especially in the 61 and 94 GHz 
(V and W) bands where the dust contribution is greatest. This model has become the standard template for 
removing the contribution of Galactic dust from microwave observations. 

To date there are no similar models that can be used to remove the polarized emission from these same 



interst ellar dust grains. WMAP polarization observations (jPage et al.ll2007t iKogut et al.ll2007t iDunklev et al 



2008bJ) at V and W-band show evidence of a dust polarization fraction ranging from ~ 1% in the direction of 
the Galactic center to a few perce nt at intermediate latitudes, consistent with the observations of Archeops 
at 353 GHz (jPonthieu et al.ll2005i ). However, due to the high noise level of the V and W-band polarization 
measurements, an estimate of the polarization angle is possible only on very large angular scales (> 4°, 
I ~ 50), and does not have high signal-to-noisc. Therefore, additional information is needed in order to 
remove polarized dust foregrounds on both large and smaller angular scales. However, this is not a limitation 
to initiating a CMBPol mission, as a dust model can be generated from the data if the mission has a rich 
enough data set in terms of frequency coverage and number of channels. This scenario occurred for both 
IRAS and WMAP, where data was used both to construct foreground templates and to measure the CMB 
spectrum. 



2.2.1 Dust Alignment 

Thermal dust polarization results from aspherical dust grains that have been aligned by interstellar magnetic 
fields. The exact grain alignment mechanism is complex so we provide only a brief description here. Basic 
alignment requires the grain's small axis, with the largest moment of inertia, to align with the spin axis, 



10 



followed by alignmen t of the the spin axis w ith the local magnetic field (e.g., Davis fc Greenstem1l951 : Whittet 



2003 



Lazarianll2003l . 120071 : iRobergd 12004). While other alignment mechanisms (e.g., mechanical alignment; 



Gold 1952) may dominate in some select environments, the above mechanism is favored in conditions prevalent 



throughout most of the ISM. Modern grain alignment theory favors grain-photon interactions, or radiative 
torques, as the mechanism by which grains achieve high-angular velocities and subsequent alignment (e.g. , 
Dolginov 1972 : Dolginoy fc Mvtrophanovlll976t iDraine fe Weingartiieilll996l Il997t lLazarian fc Hoanel 12007 ; 



Hoang fc Lazarianll2008h . The result of these mechanisms is to align the longest grain axis perpendicular to 



the local magnetic field. All grains will emit, and absorb, most efficiently along the long grain axis. Therefore, 
polarization is oberved perpendicular to the aligning field in emission, but parallel to the field in absorption 
or extinction. 

Polarized thermal emission at A = 60 - 1000 /im has been ob served in a wide variety of dusty Galactic ob- 



jects and typically has a pol a rization fraction of 0.5 - 10 % fe.g.. lDotson et alJbOOOtlCurran fc Chrvsostomou 



2007 



Matthews et al.ll2008t botson et a"l1l2009h . To date, polarization from the extinction of background 



starlight at near-visible wavelengths, believed to be due t o differential absorption, has been o bserved for more 
than 10,000 stars, typically with values ~ 1 - 4% (e.g., lHeilesll2000t iBerdvugin et al.ll2004r). Both types of 



observations have been used to study magnetic field structure throughout the Galaxy (e.g. 



200 



199 



Chuss et al 



Schleuninal Il998t iMatthews et~ZI boOlUFosalba et all l200i IZweibel fc Heilesl Il997t IPerevra fc Maealh. 



2003 



ia.es 



2007h a s well as to investigate the alignment mechanism an d physical dust characteristics (e.g. jHildebrand fc Dragovanl 



Aitken 



Fraisse et al 



ll99d|Hildebrand et alfl999ljwirittet et al.ll200ll ). For further discussion on grain alignment see 
d2008h . 



2.2.2 Polarized Dust Emission 

Current instrument sensitivity limits preclude a comprehensive characterization of the properties of polarized 
dust emission at millimeter wavelengths. Existing me asurements of diffuse regions are limited to large beams 



(intri nsic or binned) usually greater than ~ 2° (e.g.. iPonthieu et al.ll2005t iPage et al.ll2007tlTakahashi et al 



20081 ). High-sensitivity measurements of polarized dust have thus far been limited to dense regions of the ISM. 
These are generally regions in the Galactic plane (with optical extinctions exceeding Ay ~ 20 magnitudes) 
which are avoided for CMB analysis. Because dense dusty regions involve multiple populations of dust grains 
in a complex magnetic field structure, simple parameterizations describing dust behavior are likely to be 
significantly different in dense Galactic regions than in diffuse regions targeted for CMB observations, so care 
must be taken when using these measurements to predict emission characteristics in diffuse regions. However, 
useful information can still be gleaned. 

Multi-wavelength observations of dense clouds in the Galactic plane at far-infrared through millimeter 
wavelengt hs show clear evidence of a wavelength-dependen t polarization in terms o f both the fractional po - 
larization (Hil debrand et al.lll999i : IVaillancourtll2002l . 120071 ) and position angle (e.g.. ISchleuning et al.ll2000n . 
In a sample of the most massive clouds, the spectrum has a minimum at approximately 350 /im (860 GHz) 
and can vary by up to a factor of 4 in the 60 /im - 1 mm wavelength range. These changes in polarization 
amplitude and position angle are attributed to multiple populations of dust grains with different tempera- 
tures, emissivities, and polarization efficiencies, as well as differences in the aligning magnetic field direction. 
Draine fc Fraisse! (2008) demonstrate the expected variation of the polarization fraction with wavelength due 
to a mixture of silicates and carbonaceous grains. 

The structure in the 50-1000 /im (300-6000 GHz) polarizat ion spectrum of de nse clouds results from the 
existence of multiple dust temperature components at 20-60 K l|Vaillancourtil2002r ) . On the other hand, lower 
frequency emission from the diffuse IR cirrus (y < 300 GHz) is expected to be dominate d by a single cold 
dust component (T ~ 10 K), in thermal equilibrium with the interstellar radiation field (jFinkbeiner et al. 
19991 ). Therefore, although multiple polarization domains are also likely to exist along most lines of sight 



in low-extinction regions (preferred for CMB analysis) far outside the Galactic plane, this indicates that 
the polarization fraction w ill be nearly frequency-independent in diffuse clouds at A > 1 mm, v < 300 GHz 
(jffildebrand fc Kirbvlbooi . 
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Latitude 


Distance 


No. Stars (%) 


/ = P/I(%) 


E(B-V) 




Total 


4114 (75) 


1.69 


0.49 


Low Latitude 


Nearby 


1451 (26) 


0.94 


0.29 


(|6| < 10°) 


Distant 


2663 (48) 


2.09 


0.60 




Total 


1399 (25) 


0.45 


0.15 


High Latitude 


Nearby 


1315 (24) 


0.42 


0.14 


(|6| > 10°) 


Distant 


84 ( 1) 


0.89 


0.26 



Table 1: Mean stellar parameters by distance and Galactic latitude (b) from iFosalba et al.l ((2002). 
'Nearby' and 'distant' denote stars within and beyond 1 kpc, respectively. The quantities in paren- 
theses denote the fraction, in percent, of all stars in the sample. 



2.2.3 Starlight Polarization 



Background-starlight polarization is only possible in regions of fairly low extinction (Ay less than a few mag- 
nitudes for near-infrared observations) where near-visible photons have sufficient mean-free-path to traverse 
the ISM. This makes it a feasible tool for inferring the Galactic magnetic field in regions of interest for CMB 
analysis. However, the extinction places a practical limit on the most distant stars for which polarization 
can be observed. The vast majority of observed high-latitude stars are within 1 kpc of the Sun, however 
at low latitude some luminous stars as far away as 2 kpc contribute to the measured polarization signal 
(jHeileslboodlFosalba et al.ll2002h . Despite this limitation, recent analyses of such measurements suggest that 
they do contain info rmation about the uniform and random components of the magne tic field on l arge s cales 
(IFosalba et al.ll2002h . In the foreground polarization model inferred from WMAP data, I Page et al.l (|2007l ) anc 
Kogut et al.l (|2007l ) showed that the polarization angle derived from starlight polarization provides a slightly 
better fit to the data than the synchrotron polarization angle. Measurements of polarization for stars of 
different distances can also reveal the 3-D distribution of magnetic field orientations averaged along the line 
of sight. 

Tables [Hand Figure[7l show an analysis from Fosalba et al. ( 2002h using the catalog of starlight polarization 
data from iHeiled (2000). This analysis includes only stars with reliable distan ce and extinction estimates, 
and small polarization measurement uncertainties (5513 stars, 60% of the lHeilesI catalog). Low-latitude stars 
have large polarization fractions and extinctions (fd = P/I » 1.7%, E(B — V) « 0.5 mag), while high-latitude 
sources exhibit significantly lower values (fd 0.5%, E(B — V) m 0.15mag). As shown in Figure [7] there is a 
strong net alignment of starlight polarization vectors with the Galactic plane (see lower panel) as well as a clear 
alignment with the spherical shell of Loop 1 as seen from the polarization vectors in local clouds (upper panel) . 
The amplitude of th e polarization vectors are also found to vary nearly sinusoidally with Galactic longitude 



(jFosalba et alJ 120021 ). The plane-of-sky polarizations are near-minimum towards lines-of-sight tangent to 



Galactic spiral arms and near-maximum towards lines-of-sight perpendicular to spiral arms. This obse rvation 
is consistent wit h a model in which the Galactic magnetic field follows the spiral arm structure (e.g. Heilesl 



1996 



Beckl 120071 ). Further reconstruction of the three-dimensional Galactic magnetic field structure requires 



complementary observations including radio synchrotron polarizatio n, IR-submillimeter dust polarization and 



Farad ay rotational measures from distant pulsars and quasars (e.g-. lHan et alJ Han & Zhang 



2007 



Ha 



n 



2007 ). New stellar polarization data from an optical survey described in Magalhaes et al\~( 20051) i s expe cted 



in the next few years, and proposals for new observations have been submitted. See iFraisse et alJ (|2008l ) for 
further discussion. 
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Starlight Polarization (5500 Stars) 
Nearby Stars (d < 1 Kpc) 




(Fosalba et al. 2001) 



180 90 270 180 

Galactic Longitude (Deg) 



Figure 7: Starlight polarization vectors in Galactic coordinates for a sample of 5513 stars. The 
upper panel shows polarization vectors in local clouds, while the lower panel displays polarization 
averaged over many clouds in the Galactic plane. The length of the vectors is proportional to the 
polarization degree and the scale used is shown i n the lower pa n el. F or starlight the direction of 
polarization is parallel to the magnetic field. From lFosalba et al.l (|2002l ). 



2.2.4 Modeling considerations 

The expectation of changing polarization with position indicates that at a given frequency, two parameters are 
required to quantify the Q and U Stokes para meters (or P, 7) in each pixel. The fractional polarization / = P/I 
is expected to be everywhere less than ~15% ( Draine fc Fraissel2008 ). but may range from less than 1% in some 



regions to 10-15% in highly polarized regions. The dust intensity can be used to provide additional constraints 
and will be better understood with Planck observations in hand. Possible models for the polarization fraction 
include inverse variation with dust column density, or with some measure of the interstellar radiation field 
intensity. The number of parameters required to model the frequency dependence of the polarized intensity in 
a given direction is not yet established, although Sec 12.2.21 indicates that optimistically one may be sufficient 
at frequencies below 300 GHz, as the dust polarization is likely to be dominated by a single component. A 
simple power law index can be modeled with one parameter in each pixel, and could take values in the range 
1 ^ (3d 3, with (3d — 1.7 fitting current temperature and polarization data at v < 100 GHz . However, 
since the dust consists of at least two components with different indices ( Finkbeiner et al. 19991 ). additional 



parameters may be required to quantify their frequency dependence. Extrapolating from frequencies above 
~ 300 GHz is likely to lead to modeling errors, and therefore bias the estimated signal if power law behavior 
is assumed. 



2.3 Free-free emission 



Free- free (bremsstrahlung) emission is due to electron-electron scattering from warm (T e « 10 4 K) ionized 
gas in the ISM. In the optically thin regime, and at radio frequencies, the brightness temperature is given 
by T b = 8.235 x lO^aT" - 35 !/^ 1 (1 + 0.08)(EM)„ T -6 rr , where T e . is the elec t ron temperature EM is the 
Emission Measure and a is a factor close to unity ( Rvbicki fc Lightman 1979t Dickinson et al. 2003f ). The 
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electron temperature in the warm ionized medium ranges from ~ 3000 to 20000 K but is ~ 8000 K at the 
Solar Galactocentric distance. The spectrum follows a well-defined power-law with a spectral index of —2.1 
showing only a small variation with temperature. At WMAP frequencies, the effective free-free spectral 
index is ~ —2.14 for T e ~ 8000K. For modeling pu rposes, a reasonab le prediction for the intensity signal 
outside the Galactic plane is the Ha map compiled by Finkbeiner ( 20031 ) with a correction for dust extinction 
(jBennett et al.ll2003l) . or the MEM map obtained from the WMAP analysis (jGold et alj|2008h . 

Free-free emission is intrinsically unpolarized because the scattering directions are random. However, a 
secondary polarizatio n signature can occur at the edges of bright f ree-free features (i.e. HII regions) from 
Thomson scattering ( Rvbicki fc Lightman 1979 ; Keating et al. 19981) . This could cause significant polariza- 
tion (~ 10%) in the Galactic plane, particularly when observing at high angular resolution. However, at high 
Galactic latitudes, and with a relatively low resolution, we expect the residual polarization to be < 1%. There- 
fore, free-free emission is unlikely to be a major foreground contaminant for CMB polarization measurements 
in clean regions of sky. 



2.4 Additional polarized components 



There are at least two additional dust emission mechanisms that have been suggested i n the literature that 
could produce a low level of po l arized emission. A number of high Galactic la titude studies (jKogut et al. Ill996at 



de Oliveira-Costa et al 



1997 



Leitch et al 



1997 



Hildebrandt et al.l 12007). as well as studies of individual 



Galactic clouds ( Finkbeiner et al.ll2004 ; Watson et al.ll2005h . have o"bserved emission in excess of that expected 
from the three foreground components discussed above, namely synchrotron, thermal dust, and free-free 
emission. This emission has been termed anomalous for the simple reason that its provenance is not completely 
understood at this time. However, it is clear that this anomalous emissio n is highly correlated with large-scale 
maps of far infrared emission from thermally em i tting dust grains (e.g 



Dobler fc Finkbeinerl l2008at iDobler et alJ 120081 ) . Bennett et al.l (|2003f ) suggested that this correlation with 



Finkbeiner 2004; Davics et al. 2006 



the WMAP data could be e xplain ed by synchrotron emission from dusty star-forming regions. Alternatively, 
Draine fc Lazarianl ( 1998a.b,|l99^) invoked microwave emission from the dust itself to explain the excess. The 
hypothesis that this emission is dust-correlated synchrotron is disfavored by WMAP polarization maps since 
they resemble the 408 MHz synchrotron map at 23 GHz, whereas the total intensity map at that frequency 
much more closely resembles the 100 fiva dust map. Dust-correlated synchrotron would be completely ruled 
out by upcoming surveys like C-BASS at 5 GHz. Such low frequency surveys will be essential for further 
tests of these hypotheses as synchrotron and microwave dust emission have divergent spectra below the lowest 
WMAP frequency at 23 GHz. 

The correlation of anomalous emission and thermal dust emission has le d to two hypotheses in whic h 
interstellar dust grains themselves are the source of this emission (see review bv lLazarian fe Finkbeinerll2003f ). 
In the first mechanis m, small (radii ~ 0.001 /xm ), rapidly rotating grains e mit electric dipole radiation at mi- 
crowave frequencies ( Draine fc Lazarian 1998al lb; Ali-Hai'moud et al. 2008 ). In the second mechanism, larger 
(radii > 0.1 /im), thermally vibrating grains undergo flu ctuations in their magne tization, thereby emitting 
magnetic dipole radiation, also at microwave frequencies ( Draine fc Lazarian 1999h . 



T he small 'spinning-dust' g rains can be stochastically heated to temperatures approaching 100 K ([Draine fc Li 
2001[ lAli-Hai'moud et alJ 120081 ) resulting in significant mid- infrared emission (A ~ 10 - 30 /im) whereas the 
larger 'magnetic-dust' grains typically emit in thermal equilibrium in the far-infrared (A ~ 50-200 /xm). 
Observations which show stronger correlations between the excess emission and the shorter- wavelength (12 
and 25 /im) IRAS bands, as compared t o the longer- wavelength bands (60 and 100 /xm), c learly appear to 



favor spinning-dust over magnetic-dust (|de Oliveira-Costa et al 



2002 



Casassus et al 



200 



6). 



Furthermore , 



correlations of the microwave excess with Ha emission peaking at ~40 GHz (|Dobler fc Finkbeinerl 12008b 
Dobler et al. 20081) also favor spinning-dust in the warm ionized medium (WIM) , as both Ha and WIM spin 
ning dust are expected to be proportional to the density squared of ions. Despite the strong evidence for 
spinning dust, emission from magnetic-dust must exist at some level, as large grains are known to exist from 
observed emission in the far infrared, and contain ferromagnetic material which is depleted from the gas- 
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phase (e.g. Savage fc Sembach 1996t Sembach fc Savage 1996 ). This is especially important for polarization 
observations as magnetic dust is predicted to be much better aligned than spinning dust. 

The small grains responsible for spinn ing-dust emission are the same g rains responsible for starlight 
extinction at ultraviolet wavelengths (e.g.. iMathis. Rumpl. fc Nordsiecklll977l ). Theoretical models of grain 
alignment predi ct spinning-dust polarizat ion could be as high as 7% at 2 GHz, falling to < 0.5% at frequencies 
above 30 GHz ( Lazarian fc Drains 2000). Observationally, the relati vely low polar i zation observ ed at UV 
wavelengths suggests that these small grains are inefficiently aligned (jWhittetl l2004t iMartinl 120071 ) and will, 
therefore, produce little polarized emission at any frequency. The larger vibrating magnetic-dust grains are 
also responsible for the thermal emission seen at far-infrared and submillimeter wavelengths. There is clear 
evidence that these grains are well aligned since their polarized thermal emission is observed to be as high 
as 10% f ^2.2.1[) . Theoretical models sugge st that magnetic dipole emission from these grains can reach levels 
as high as 40% (jDraine fc Lazarian! 1 19991 ) . Additionally, the polarization angle is predicted to exhibit 90° 
flips as a function of frequency within the ~ 1 - 100 GHz range; at high frequencies the position angle is 
perpendicular to the aligning magnetic field (parallel to the angles of polarized spinning-dust, thermal dust, 
and synchrotron emission), but parallel to the field at lower frequencies. 

Regardless of the exact emission mechanism, the polarization will need to be explored and accounted 
for in foreground cleaning, especially since the morphology and frequency dependence of the emission is 
not par ticularly well known. The 23 GHz WMAP map, as well as targeted observations o f specific Galactic 
clouds (jMason. Robishaw. fc Finkbeinerll2008l : iBattistelli et alJl2006t iDickinson et al.1120071 ) suggest that the 
polari zation of the anomalous emission is low, £ 5%, with a measurement of 3.4j^'g% by IBattistelli et al. 
(2006). This is consistent with the models of polarized spinning-dust discussed above, but cannot rule out 



significant contributions from polarized magnetic-dust emission. It is worth noting however that the challenges 
of addressing these anomalous mechanisms are not insurmountable, and similar challenges have been addressed 
successfully with previous experiments including FIRS, DMR, FIRAS, and WMAP. There are many ways 
to approach and remove foregrounds provided that low noise observations are made at a sufficient range of 
frequency channels at 20 - 1000 GHz. Prior to the launch of a CMBPol mission, ground based surveys will 
cover the lower frequency range. In addition, complimentary targeted searches of specific Galactic clouds will 
continue to constrain the polarization properties. 

Finally, it is possible that in addition to these dust emission mechanisms, there may also be other exciting 
discoveries awaiting us that at present are not in our models. 



2.5 Simulated foreground maps and spectra 

By drawing together current observations and theoretical predictions, we generate simulated maps of the 
polarized Galactic emission with pixels of side 4 and 0.5 deg reefl In anticipat ion of the Planck satellite 
mission, the 'Planck Sky Model' (PSM) has been developed (jLeach et al.l 120081 ). drawing on observations 
including WMAP, IRAS, and 408 MHz radio data. The Planck Working Group 2 have developed this model 
and allowed us to use components of it for this work. For initial forecasting purposes we generate a set of 
foreground maps with synchrotron and dust components, each with power-law spectral indices. This physical 
behaviour is over-simplified but useful for preliminary tests. The maps with 4 degree pixels are smoothed 
with a 7 degree Gaussian beam, and those with 0.5 degree pixels with a 1.5 degree beam. 
The synchrotron component is formed at each frequency v using: 



Q sy „ci>, n) - (v/v s y^Q s (n), U synch (v, h) = {v/v s f'^U s {n) 



(4) 



For initial studies, the maps for Q s 
(20080. They are formed using 



U s and (3 S are generated from the model of iMiville-Deschenes et al 



Qs{n) = fs9s(n) cos(27 s (n))J s (n) U s {h) = f s g s (h) sm(2"/ s (n))I s (h), 



(5) 



use HEALPix iV sidc =16 and 128. 
2 As implemented in the PSM vl.6.2. 
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Synchrotron P (90 GHz) Dust P (90 GHz) 




Figure 8: Simulated maps for synchrotron (left) and dust (right) polarization amplitude in antenna 
temperature, polarization angle, and spectral indices at 90 GHz, drawn from the Planck Sky Model 
and described in the text. 



where I s is the pred icted synchrotron intensity at v s = 33 GHz extrapolated from observations at 408 MHz 
(jHaslam et al.lll981l ). The geometric suppression factor g s , the polarization angles j s , and the polarization 
fracti on f s are obtaine d using a magnet ic field model constrained usin g the WMAP 23 GHz polarization 
data 1 Page et al. 2007 ) and described in Miville-Deschenes et al. ( 2008h . The geometric suppression factor 
accounts for spatially varying depolarization due to line-of-sight integration through the Galactic magnetic 
field, resulting in a lower observed polarization fraction than at the point of emission. The factor does not 
account for beam dilution. 

The simulated spectral index map j3 s is obtained by first estimating the synchrotron intensity at 23 GHz, 
using the polarization maps and the same Galactic magnetic field model described above, assuming an emitted 
polarization fraction of f a = 0.75. A power law is then fit to the 408 MHz Haslam intensity data and this 
estimated 23 GHz intensity map. This ignores any spectral break that one might expect in that frequency 
range, but a scale dependent spectral index can be added to the simulations. The polarization amplitude, 
angle, and power law spectral index are shown in Figure [H 

A simple parameterization for the synchrotron emission is 



t(t + i)cf nch /27r = A s (v/v ) 2 ^(£/e y 



(6) 



although this oversimplifies the true behavior. Estimates for these quantities are A s ~ 0.012 iiK 2 , 0. ~ —3.0, 
and m s ~ —0.6 at v s = 90 GHz and Iq = 10, based on WMAP observations by I Page et al.l (|2007h . Figure 
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Figure 9: Compilation of measured and predicted E-mode and B-mode amplitudes at 90 GHz for 
synchrotron and dust emission compared to the lensed scalar (long dash blue) and tensor (solid blue) 
CMB signal with r = 0.085 and r = 0.01. Spectra are shown in antenna temperature units for a sky 
fraction of 75%, and 3%, in the estimated lowest foreground region. The synchrotron contribution 
at large scales is constrained by WMAP observations. The dust contribution, and the behavior of 
both components at £ ^ 30, is uncertain. The red lines show two realistic levels for the predicted 
dust spectrum, with 5% and 1.5% average polarization. Only the lower level is shown for the 3% 
sky fraction. Recent simulated dust maps by the Planck collaboration give a B-mode spectrum that 
lies within these limits, and has a flatter spectrum. 



[9] shows this spectrum as a function of angular scale, compared to the primordial CMB signal with tensor- 
to-scalar ratio r = 0.01. The estimated synchrotron amplitude (given by the square root of the spectrum) is 
about three times bigger than the CMB signal for r — 0.01 at I = 100 over 75% of the sky at 90 GHz, and 
about ten times bigger than the CMB signal at I = 4. Over the cleanest 3% of the sky the signal is expected 
to be comparable to the CMB signal for r = 0.01 . 

There is far more uncertainty in the polarization of the thermal dust component, with about an order of 
magnitude uncertainty in the amplitude, ranging from ~1% to 10% polarized. To make simulations, the dust 
component is formed in a similar way to the synchrotron, using 

Qdu S x{^n) = {u/u d f^Q d {n), U dust (u, n) = {v/v d f^U d {n). (7) 

The amplitudes and spectral indices are taken from the Planck Sky Model, with the amplitudes formed using 

Qd{n) = fd9d(n) cos(2j d (h))I d (n), U d (h) = f d g d (n) sm(2-f d (n))I d (h). (8) 



Here I d is the dust intensity at 94 GHz extrapolated from IRAS by iFinkbeiner et alj ((1999) . In the maps 
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we adopt as fiducial estimates, the polarization fraction fd is set to be 5%, and the polarization angles 
and geometric suppression factor are calculated using the Galactic field model described for the synchrotron 
simulations. This corresponds to the Planck Sky Model vl.6.2, and assumes that the dust and synchrotron 
have the same po larization angle. This is not correct along certain lines of sight in the Galaxy (see e.g., 
Page et al.l (|2007l )). but is used as a first approximation. The resulting map is only 1-2% polarized due to 
depolarization effects. In more recent investigations for Planck, a higher polarization fraction of 12% at 
emission is assumed, resulting in an observed polarization fraction of ~ 4%. For the frequency dependence 
we use power law indices calculated using FDS Model 8 between 84 and 94 GHz, with an index j3d = 1.65. 
Figure [5] shows the polarization amplitude, angle, and index for these maps at v& = 90 GHz. A more realistic 
model would have multiple polarized dust components which can be modeled by adding their contributions 
to Q and U given in Eqn. [8] A two component model is used in the current Planck Sky Model. 

Figure [9] shows the angular power spectrum at 90 GHz for the maps outside the P06 sky cut, fit with a 
power law that flattens at large scales. This corresponds to ~ 1.5% polarization fraction. The power spectrum 
currently estimated by the Planck collaboration (PSM vl.6.4), has a higher polarization (3-4%) and a flatter 
scale dependence. This mode l account for more of the turbulent effects in the Galactic magnetic field that 
tend to flatten the spectrum ( Prunet et al. 19981 Prunet fc Lazarian 19991 ). For comparison we also plot 5% 
of the FDS dust intensity. Parameterizing the dust as 



£{£ + l)Cf ust /27T = A d (v/v Q ) 2 ^(e/£ y 



0) 



fails to correctly describe the large scale behavior, but is a useful starting point. Estimates for these quantities 
over 75% of the sky fall within typical ranges 1.5 < (3d < 2.5, -0.5 < m < 0.5, and A d < 0.8 fiK 2 for £ = 10 
and Vq — 90 GHz. The simulated maps with 1-2% polarization fraction have Ad = 0.004 /iK 2 at £ = 10, 
(3d = 1-65, and typically m — —0.5. For comparison, at 90 GHz, 5% of the FDS in t ensity corresponds to 



A d ~ 0.1 /iK 2 at £ — 10. An upper limit of 15% is suggested by Draine fc Fraissei |2008h . For a smaller 



patch of sky, taking a circle of radius 20 degrees in a region of low foreground intensity, the predicted B-mode 
amplitude is significantly lower, with A d = 6 x 10~ 5 /iK 2 , and A s ~ 0.0002 /iK 2 at £ — 100 and v = 90 GHz 
for the simulated maps. 

Figure [10] shows the total predicted power in the diffuse foregrounds maps as a function of frequency for 
£ — 80— 120j. This is based on the PSM model (vl.6.3), with a dust polarization fraction of onl y 1-2% percent. 



The C MB amplitude is for a standard ACDM model as derived from WMAP five-year results (jDunklev et al 
2008a[) . including tensor fluctuations with r = 0.01. Extragalactic sources are not included. The maximum 
ratio of CMB to foregrounds occurs at ~ 100 GHz and moves to slightly higher frequencies for cleaner patches 
of sky. At 100 GHz and £ = 100 the foreground signal is ~ 5 times larger than the B-mode signal for r = 0.01 
outside the Galactic plane (in the WMAP Kp2 mask), but in a clean patch of radius 10 degrees the foreground 
signal is below the CMB signal. However, there is still about an order of magnitude uncertainty in the total 
polarized dust contribution which should be resolved with observations from Planck, and from ground and 
balloon-based experiments. 



3 Methods for polarized CMB estimation 

In this section we briefly review and contrast methods commonly used to estimate the polarized CMB signal 
from raw sky maps. We focus on those that have been applied to current polarization data, and that can 
be used to forecas t the performance of a future m ission. For a description of additional methods see e.g., 
Leach et al.l |2008l ) ; IPeTabrouille fc Cardoso! |2007h . 



The power spectra are calculated with PolSpice v2.5.7 (see e.g. IChon et alj (|200J)). 
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Figure 10: Predicted diffuse foreground power at angular scales t = 80 — 120, as a function of 
frequency and sky coverage, compared to a CMB signal with r = 0.01. From top to bottom are the 
TT, EE and BB power spectra in antenna temperature. The CMB is constant in thermodynamic 
temperature and thus decreases with frequency in these units. Left: Total power for different sky 
coverage: full-sky, \b\ > 10°, \b\ > 30°, \b\ > 50°, and a clean circular patch of radius 10° centered on 
(I, b) = (240°, —70°). Right: Ratio of the total diffuse foreground power to the CMB. The maximum 
ratio occurs at ~ 100 GHz and moves to higher frequency for cleaner patches of sky. 
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3.1 Template cleaning 

A simple foreground cleaning technique assumes that a map of the microwave sky T(p,v), at some pixel p 
and frequency is, can be described as a superposition of fixed spatial templates X(p) and noise n(p, v), 



T(p,v)=n(p,v) + J2<Xi(») x i(p), (10) 

i 

where the coefficients oiiiv) describe emission traced by the i th spatial template at observing frequency v. 
Provided that the noise covariance matrix (which typically includes both instrument noise and the CMB) is 
known, Eqn. [10] can be solved for the template coefficients a%{y). The noise in the foreground-cleaned map, 

T c (p,v)=T{p,v)-Y,(*i{v)Xi(p) (11) 

i 

is nearly unaffected by template cleaning and retains the properties of the noise in the original map. This is an 
important consideration for applications requiring accurate propagation of noise into, for instance, estimates 
of cosmological parameters. The amplitude of the noise in the cleaned map increases as more parameters are 
fit; however, for most applications there are many more pixels than fitted parameters so the resulting increase 
is negligible. A typical sky map with 10 5 pixels could simultaneously fit 1000 spatial templates and suffer 
only a 1% increase in the noise amplitude after cleaning. 

Template cleaning has a number of other attractive features. Since each template is fit to all map pixels 
simultaneously, the technique makes full use of the spatial information in the template map Xi(p). This is 
important for the non-stationary, highly non-Gaussian emission distribution typical of Galactic foregrounds, 
allowing cleaning of the entire map even if the foreground emission is faint compared to the noise in each 
individual pixel. The method also allows multiple template maps to be fit to a single frequency channel, as 
opposed to pixel-by-pixel techniques which generally require at least one frequency channel per foreground 
component to be fit. In addition, template cleaning is insensitive to spatial correlations between various 
foreground components. Correlations between the template maps will create non-zero covariance in the fitted 
parameters o^; however, this simply complicates the identification of the associated emission Tj = ctiXi with 
the emission mechanism traced by that template, but does not bias the total foreground estimate summed 
over all templates. 

Several cautions apply to template cleaning. The method assumes that the spatial and frequency depen- 
dence of each emission component can be described using separable functions for the spatial and frequency 
dependence, T(p, v) = X(p)f(i>). Many physical emission processes violate this assumption at least to some 
extent. The spectral index of synchrotron emission, for instance, is known to vary with position on the sky, 
while dust is known to require at least two temperature components whose ratio also varies across the sky. 
Over sufficiently large separations in frequency, such spatial variation in the frequency scaling can significantly 
alter the spatial distribution of the affected emission component. In addition, the technique requires a separate 
template for each spatially distinct foreground and is thus ineffective as a 'blind' test for foregrounds whose 
spatial distribution is not well matched to one or more templates. Template cleaning works best when the 
template map is signal dominated. Noise in the template map will be aliased into the cleaned map and could 
become significant in clean regions of the sky if the template noise is sufficiently large. While this effect can 
be estimated using Monte Carlo techniques, some care in template selection and sky cuts should be exercised. 

Figure QT] shows a selection of templat e maps widely use d for CMB analyses. Synchrotron emission is often 
fit using the unpolarized 408 MHz survey ( Haslam et al.|[l981 ). Due to concerns over possible spatial variations 



in the synchrotron spectral index, the WMAP team has also constructed a sy nchrotron-dominated template 



using the difference of the two lowest frequency channels (jHinshaw et al.ll2007| ). The CMB anisotropy cancels 
exactly in such a difference map, leaving a dominant synchrotron signal with smaller contributions from free- 
free and dust, as well as instrumental noise. Provided the template fit (Eqn. I10|) includes separate templates 
for free-free and dust, the resulting linear combination of templates will remove all three foregrounds. As 
described in Section [2T3l free- free emission from the warm ionized interstellar medium may be traced using Ha 
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Figure 11: A sample of template maps commonly used to clean polarized and unpolarized maps. 



emiss ion from the same ionized gas, subject to non-negligible corrections for optical extinction by interstellar 
dust ( Finkbeiner 20031 : Schlegel et al. 1998h . Dust emiss ion is commonly traced using a model of millimeter- 
wave dust emission based on the COBE and IRAS data ([Finkbeiner et al.lll999l ). described in Section |2~21 and 
shown in Figure [6l 

Templates for polarized emission are harder to come by. Free-free emission is largel y unpolarized, so 
the WMAP 23 GHz polarization map can be used as a polarized synchrotron template (jPage et al.l 12007 : 
Gold et al.ll2008l ). No high signal-to- noise ra tio tracer of polarized du st emission currently exists. A polarized 
dust template may be constructed using the iFinkbeiner et al.l (|1999i ) unpolarized dust model convolved with 
optical measurements o f the dust polari zation angle and a geometric model of the Galactic magnetic field 



along each line of sight (Pa ge et al 



2007). 



Template cleaning is widely u sed for CMB an alysis. As an example, Table [2] shows the coefficients 

2008h . Separate cleaning is performed for the polarized and 



fitted to the WMAP 5-year data (jGold et al 



unpolarized maps. The unpolarized maps are fit o yer the high- l atitud e sky using three spatial templates: the 
syn chrotron-dominat e d K-K a difference map, the IFinkbeiner! ([2003) Ha map corrected for extinction, and 
the IFinkbeiner et alJ (|1999l ) model 8 dust map evaluated at 94 GHz. Since the "synchrotron" template uses 
both the K and Ka channels at 23 and 33 GHz, the template cleaning is applied only to the 3 high-frequency 
unpolarized bands. The polarized data, at reduced pixel resolution, are fit using the WMAP 23 GHz map as 
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Freq 
(GHz) 


Unpolarized Coefficient 
Synchrotron Free-Free Dust 


Polarized Coefficient 
Synchrotron Dust 


33 








0.317 0.017 


41 


0.24 


1.00 


0.20 


0.177 0.015 


61 


0.06 


0.65 


0.47 


0.060 0.037 


94 


0.00 


0.40 


1.26 


0.045 0.082 



Table 2: WMAP Template Coefficients an(u) 



a synchrotron tracer, plus the iPage et al.l (|2007l ) template for polarized dust. 

Template cleaning is remarkably efficient at compressing complicated Galactic foregrounds into a handful 
of parameters. Nine parameters describe the unpolarized foreground emission from some seven million input 
pixels, while 8 parameters describe the polarized emission from nearly 20,000 pixels. This simple cleaning 
technique produces results in agreement with considerably more complicated pixel-by-pixcl methods which 
include additional astrophysical information. 

Template cleaning can also lead to surprising results. Figure IT21 shows the template coefficient a(u) for 
a template map of thermal dust emission at 240 fim fitted to the COBE/DMR microwave maps at 31, 53, 
and 90 GHz fcoeut et alJll996al lbh. At sub-mm wavelengths the coefficients follow the expected behavior for 
thermal dust emission, but show an unexpected increase in emission at mm wavelengths. This "anomalous 
emission", correlated with the spatial distribution of thermal dust but with spectral index (3 ~ —2.2 (in 
antenna temperature) instead of the values 1.6 < (3d < 2 expected for thermal dust at mm wavelengths, has 
since been confirmed in a number of data sets. This is the same emission that is discussed in Section [2.41 and 
although its physical origin is still unclear, it is efficiently removed by template cleaning. 

Template cleaning for CMB analysis to date has operated in two distinct regimes. Unpolarized fore- 
grounds at high latitude are typically fainter than the CMB anisotropy, requiring only a light cleaning of 
weak foregrounds. Polarized foregrounds are comparable to or brighter than the CMB polarization over most 
of the sky. Current detections of E-mode polarization, at amplitudes comparable to the foregrounds, thus re- 
quire modest cleaning of bright foregrounds. A search for primordial B-mode polarization will push template 
cleaning into a third regime, requiring deep cleaning of bright foregrounds. Detecting a signal with tensor to 
scalar ratio r = 0.01 at multipoles £ < 10 (angular scales 9 > 20°) requires reducing the foreground emission 
by a factor of roughly 20. 

Happily, recent history suggests that such deep cleaning is feasible and should be possible with a future 
mission. The cosmic infrared background (CIB) observed at wavelengths 125 < A < 2000 /im results from the 
integrated infrared emission of galaxies at redshift z ~ 3. Galactic emission at these wavelengths is brighter 
than the CIB and has a similar spectrum. Template cleaning of the COBE/FIRAS data using only 3 spatial 
templates (the 158 /im line from ionized carbon, the 21 cm line from neutral hydrogen, and the square of 
the 21 cm line intensity) successfully reduced a bright, complicated f oreground by a fact or of 10 to enable a 
statistically significant detection of the CIB intensity and spectrum (jFixsen et al. 1998 ). Furthermore, the 
ionized carbon template was derived from the FIRAS data itself, providing an example of a rich data set being 
used for internal cleaning. Detecting B-mode polarization in the CMB would require only modestly deeper 
cleaning, and would have the distinct advantage that the CMB frequency spectrum is known a priori. 

Template cleaning has been demonstrated to reduce foreground emission by at least a factor of 10. The 
ultimate limit for this technique is not known. Given the large number of existing maps for various Galactic 
emission sources and the small noise penalty to be paid for fitting multiple templates, an interesting scenario 
would be to attempt a simultaneous fit of polarization data to as many templates as possible, including (as 
was done for the CIB) higher powers of individual templates. There will be additional templates available 
in the next five years, at low frequencies (5-20 GHz) from the C-Band All-Sky Survey (C-BASS), GEM-P 
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Figure 12: The frequency dependence of emission correlated with thermal dust shows an anomalous 
increase at centimeter wavelengths. Despite uncertainty in the phys ical origin of this c orrelated 
emission, it may efficiently be removed using template cleaning. From lKogut et al.l (|1996al ). 



( Barbosa et al. 200(jh . and COFE ( Leonardi et al. 2007 ). and at higher frequencies from the Planck satellite 



3.2 Parametric methods 

Bayesian parameter estimation methods have begun to be applied to the problem of estimating the CMB 
signal, given observations of the sky at multiple frequencies. In contrast to template cleaning, which produces 
a cleaned map at each frequency, these methods use the multi- frequency information to estimate a single CMB 
map. The methods exploit our knowledge of the frequency dependence of different components. The signal 
is parameterized by a model, expressed as the sum of CMB and foreground components at each frequency. 
Current methods are split into two schemes. In the first approach the pixel Q and U CMB Stokes parameters 
are estimated in each pixel. In the second approach the CMB angular power spectrum is simultaneously 
estimated with the map parameters. In this framework, templates comprising the 'best-guess' for the emission 
of a particular foreground component can be included as priors. One of the main benefits of this method is 
that errors due to foreground uncertainty are rigorously propagated into the estimated CMB maps and power 
spectra. One of the drawbacks is that polarized components that are not in the parameterized model, or that 
are incorrectly described by the model, will not be properly accounted for. These methods are also typically 
computationally demanding. 



3.2.1 Estimation of CMB maps 

Eriksen et al] |2006h . bold et all |2008h . bunklev et all |2008bl ) and Istompor et all |2008h describe a CMB 
component separation approach based on standard Bayesian parameter estimation techniques. A paramet- 
ric model is used for each signal component (including CMB, synchrotron, and dust) and the probability 
distribution for the parameters is estimated. 
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Using Bayes' theorem, the posterior probability is given by 



P(9\d) oc P(d\6)P(9) = £{0)P(0), 



(12) 



where C(9) = P(d\9) is the likelihood of the observed maps, and P(9) is a prior for model parameters 9. The 
likelihood is given by 

-21n£ = £[d„ - S v (e)] T N- 1 [d w - S„(0)], (13) 
i/ 

where N„ is the noise covariance at each channel v, and the data contains some combination of I, Q, and U 
in map space. The model, S v {&), is minimally parameterized as the sum of CMB, synchrotron emission, and 
thermal dust emission components: 



= ol c , v A c + ot SiV A a + a dyV A d . 



(14) 



For the CMB the coefficients are already known: the emission is black-body. For dust and synchrotron 
emission the coefficients are often constrained to give power law spectral indices, with a v {n) = (v / 'v§)^ n ' . 
The amplitudes maps Aj arc defined at pivot frequencies vq. The goal is to estimate the CMB map A c and 
its associated covariance. 



The joint distribution for the amplitudes and s pectral indices cannot be sampled directly, so i n Erikse n et al 



2002 



(2006 ) the MCMC Metrop olis-Hastings algorithm |Metropolis et al.ll953Uchristensen et alj200ll:lLewis fc Bridle! 



Dunklev et al.ll2005l ) is used to map out the the posterior probability distribution for the parameters 9, 



using intensity data (J). It has also been applied to polarization data (using /, Q, U), in studies for the Planck 
satellite. In this case the noise is assumed to be diagonal, so each pixel is fit individually. The full parameter 
set, including component amplitudes and spectral indices, is initially fit pixel by pixel in low-resolution maps. 
This provides both best-fit values for each parameter and its associated uncertainty. The second sampling 
step is performed at higher resolution. All non-linear parameters are set at the best-fit values obtained in the 
low-resolution sampling, the parameter fields are smoothed spatially, and high-resolution amplitude maps for 
each component are then obtained analytically. The residual uncertainty in the measurement of the CMB 
emission in each pixel after foreground signals have been removed is estimated using the multiple bands. The 
technique is particularly appropriate given that Galactic emission anisotropy power dominates on large scales. 
The code 'FGFit' provides an impleme ntation of this met hod and is also a useful forecasting tool. A similar 
method has also been implemented by Gold et al. ( 20081 ) and applied to the 5-yr WMAP temperature and 
polarization data. 



An alternative method to estimate polarized CMB maps is described in iDunklev et alJ (|2008bl ). The 
principle is the same, but the sampling method is modified, with a Metropolis-within-Gibbs MCMC method 
used to map out the joint distribution p(9\d) — p(A,a\d). The Gibbs sampler is a special case of the 
Metropolis-Hastings MCMC algorithm. Specifically, suppose the target distribution of interest is P(A,B), 
and we want to generate a joint sample (A, B) from this distribution. This can be achieved by alternately 
sampling from the two conditional distributions, 



A i+i 
B i+i 



P{A\B l ) 
P{B\A l+1 ) 



(15) 
(16) 



Here, the left arrow indicates sampling from the distribution on the right hand side. In this method the param- 
eter distribution is sliced into a conditional distribution for the amplitude maps, p(A|/3, d) for fixed non-linear 
parameters, and a distribution for the spectral indices p(/3|A, d), for fixed amplitude. The first distribution 
is a Gaussian and can be sampled directly. The second distribution is sampled using the Metropolis-Hastings 
algorithm. An advantage of this approach is that the noise matrix does not have to be assumed diagonal, and 
spectral indices can for example simultaneously be estimated in larger pixels than the amplitudes. P r iors are 
imposed on the foreground emission. For example, for application to WMAP data, bunklev et~aH |2008bl) 
assume that the spectral indices in Q and U are the same in a given pixel, and that in a map with 768 pixels, 
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the spectral indices are only allowed to vary in iVj = 48 larger pixels. In regions of low signal-to-noise priors are 
imposed on the synchrotron and dust spectral indices. A typical choice is a Gaussian prior of (3s — —3.0 ±0.3 
on the synchrotron emission index, and (3d — 2 ± 0.3 (or broader) for the dust emission index. 

While the methodology is different for these two methods, the principle is the same, and the resulting 
product is a CMB map and covariance matrix estimated from the marginalized posterior distribution. The 
likelihood of the estimat ed maps, given a theoretical angular power spectrum, can then be computed using 



the method described in I Page et al.l (|2007l) . By varying only the tensor-to-scalar ratio r, and calculating the 
likelihood at each value of r, this method can be used to estimate limits on the tensor-to-scalar ratio including 
foreground uncertainty. 

3.2.2 Estimation of CMB power spectra 

Rather than estimating the marginalized maps, and then estimating the CMB power spectrum, it is possible 
to simultaneously estimate both. Here the main goal is the joint posterior distribution P(s, CV, #f g |d), where 
s is the CMB sky map (denoted A c in the previous sub-section), Ci is the CMB power spectrum, 8{ g denotes 
the set of foreground parameters, and d represents the data. From this joint distribution, one may extract any 
number of marginal distributions to obtain estimates for each parameter individually, of which perhaps the 
most important is the marginal power spectrum distribution P(Ci |d). From this one may extract cosmological 
parameters using a standard MCMC code such as CosmoMC. 

The joint distribution P(s, Cg, 0f g |d) involves several millions of correlated parameters, and this poses 
a serious computational challenge. Brute force evaluation is impossible. However, in recent years a new 
Monte Carlo method based on the Gibbs sampling algorithm has been developed for precisely thi s purpose 
(jjewell et al.ll2004 IWandelt et~al]|2004 lEriksen et alll2004 lLarson et alll2007t lEriksen et aDl2007bh , and this 



allows the user to draw samples from the full joint distributi on in a computational l y efficient manner. Ap pli 



cations to real- world data include the analyses presented by lO'Dwver et al.l (|2004l ): lEriksen et al. 

In order to use the Gibbs sampler to sample from the joint CMB posterior, one must simply be able to 
sample from all corresponding conditional distributions, 

cj+i ^ p{c,\s\e\ g ,d) (17) 

s' 1+1 ^P(s|q +1 ,^ g ,d) (18) 
^ P(9 fg \Cl +1 , s*+\d). (19) 

Then, after some burn-in period, the set s 1 , 6\^} will be drawn from the desired joint distribution. Fortu- 
nately, sampling from each of these conditionals is straightforward, as the distribution for Ce in Eqn. [T7]is an 
inverse Wishart distribution, the distribution for s in Eqn. [18] is a multivariate Gaussian, and the distribution 
for 9{ g is in general given by a x 2 . All of these have well-known sampling algorithms (see, e.g., Eriksen et al. 
2008a for a full review). 

As in Section l3.2.H the foreground model may be specified quite freely. For application to current data, the 
focus has been mainly on models with a small set of parametric foreground spectra with individual parameters 
for each pixel. Typical examples include power-law spectra (S u — Av 13 ; A and (3 free), one- or two-component 
thermal dust spectra, SZ spectra and arbitrary tabulated spectra. Additionally, Gaussian and/or uniform 
priors may be specified for each pixel. In the future, complementary sets of parameters describing a spatial 
power spectrum for each foreground component may be used, a feature which will be useful for probing the 
low signal-to-noise regime. 

A major advantage of the Gibbs sampling approach is, because of its conditional nature, its ability to 
integrate additional sources of uncertainty without modifying the existing framework. Three specific examples 
would be gain calibration uncertainties, beam uncertainties and noise estimation errors. By adding another 
conditional sampling step for each of these to the above Gibbs chain, the corresponding uncertainties would 
be seamlessly propagated through to all other parameters, without requiring any modifications to the main 
sampling scheme. 
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3.3 Blind component separation 



In both the template and parametric methods, the foreground model is the sum of specific astrophysical 
components, and its construction relies on knowledge of their spatial behavior and/or frequency dependence. 
The goal of blind or semi-blind component separation methods is to separate the observed sky maps into 
a CMB component and a set of foreground components, with few assumptions about the physical behavior 
of those foreground components. Here we describe two methods, Internal Linear Combination (ILC) and 
Independent Component Analysis (ICA). 



3.3.1 Internal Linear Combination 



The ILC technique in spherical harmonic 



CMB anisotropies in lTegmark fc Efstathiou 



five-year temperature data IIHinshaw et al 



Hinshaw ct al. 2007; Gold et al 



s pace was first described as a way to remove foregrounds from 
It was applied successfully to the WMAP one, three, and 



spherical harmonics (jAmblard et al 




Tegmark et "all l2003t Ide Oliveira-Costa fc Tegmarkl 12006 



2008) , but it can also be applied to polarized data, using the E or B mode 



2007). The CMB signal is estimated as a linear combination of the 



observed a\ m in each frequency band i, 



frcq— ? 



w\a l7) 



(20) 



with weights Wi chosen to minimize foreground contamination. For polarized CMB observations, the signal 
at each frequency can be decomposed as 



Clr, 



Ini 



(21) 



where c, s, d, and n stand for the CMB, synchrotron, dust, and noise. The number of foregrounds can be 
increased, provided that the number of channels is also raised. To choose the weights that best extract the 
CMB signal, the power spectrum 

(Km | 2 } = w/C/Wf, (22) 

is minimized to optimally remove the foregrounds, where C 1 / = ((d| TO ) a £ m )- The weights are constrained by 
wj - e = 1, where e is a column vector of all ones with len gth equal to the numbe r of channels. This condition 
ensures that the CMB signal is retained. As derived in begmark etaD fcooah . the minimized weights are 
given by 

J- • (23) 



w f = 



3"^ C n 



The estimated spherical harmonics are then used for cosmological analysis. The advantage of the ILC al- 
gorithm is that there are no assumptions made about the frequency dependence of the foregrounds. The 
algorithm maximizes the signal to noise of the recovered CMB, and accounts for differences in beam size 
between frequencies. The algorithm is computationally very fast, allowing multiple Monte-Carlo simulations 
to be performed to propagate systematics and compute errors. A disadvantage of the method is that the 
foregrounds are assumed to have a spatially invariant frequency behavior, which is not physically realistic. 
This can be relaxed by splitting the full sky maps into several areas. Further drawbacks are that the total 
power spectrum is minimized rather than just the foreground power, and it is difficult to rigorously propagate 
errors into the estimated CMB maps. 



3.3.2 Independent Component Analysis 

As with the ILC method, ICA components are linear combinations of the input sky maps, but the weights 
are chosen using an alternative method. Two examples are 'altICA' and Spectral Matching ICA (SMICA). 
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The altICA method is based on an algorithm by iHvvarinenl (|l999| ). and is described in iMaino et al 



(20020. The weights given to the sky maps for each component are estimated by maximizing the non- 



Gaussianity of the component maps. This idea is based on the central limit theorem, which states that a 
variable that is a mixture of independent variables is more Gaussian than the original ones. By transforming 
the observed sky maps into a set of components such that the Gaussianity of the variables is reduced, the 
resulting maps should be more independent. The algorithm has been implemented as a CMB cleaning pro- 
cedure, because the CMB and diffuse foregrounds are expected to be statistically independent. However, the 
non-Gaussianity of a given foreground component may vary spatially, so this decomposition is not necessarily 
astrophysically motivated. 

In this method, the degree of non Gaussianity is quantified by the neg-entropy statistic, which is more 
robust than many other statistics in the presence of noise. The neg-entropy for a map y (containing /, Q or 
U, for a given set of weights), is defined as 

neg-entropy(y) = H(y G ) - H{y) , (24) 

where H(y) = — J p(y) logp(y)dy is the entropy associated with the distribution p(y), and y^ is a Gaussian 
variable with the same covariance matrix as y. This method can be extended to include a priori knowledge of 
the frequency scaling of one or more components. The method has been applied to real and simulated inten- 



sity data (including BEAST, COBE and WMAP) and simulated polarization data (see IMaino et all (|2007j ) 
and references therein) . The performance relies on the assumption of statistical independence for CMB and 
foregrounds, and operates most successfully with high resolution data. 



In the SMICA method (jDelabrouille et al.l 12003 : ICardoso et all 120081 ). the CMB and foreground power 



spectra are estimated using a parametric model o f the angular power spectra of the components. For appli- 
cation to Planck simulations ( Cardoso et al. 2008h . the angular spectra and cross-spectra, C 1 / , measured at 



N v frequencies, are parameterized as the sum of CMB, foregrounds, and noise: 

Cf = C £ MB >v + pf + N*i , (25) 

with C 1 / — l/(2£+ 1) X) m ( a im)^ a im f° r frequencies i and j. In practice the spectra are binned in N q bins 
in £, to give spectra C q J . The foreground model is parameterized by D possibly correlated angular power 
spectrum templates, with 

F q = A T P q A. (26) 

Here A is a matrix of size D x N u and P q is a symmetric matrix of size D x D. This gives a total of 
N v x D + N q x D(D + l)/2 foreground parameters. SMICA allows foregrounds to be correlated, which 
is useful for modeling Galactic components that are naturally correlated by the geometry of the Galactic 
plane and the effect of the Galactic magnetic field. To capture the complexity of the Galactic emission in 
polarization, at least two templates are needed for synchrotron and dust emission. In practice, the number of 
components required to capture the physical behavior depends on the frequency range and noise levels of the 
experiment. To study CMBPol missions, D = 4 templates are used for the Galactic emission, which captures 
synchrotron and dust emission, and corrections due to variable emission laws. This parameterization of the 
foregrounds is considered 'semi-blind' as it allows only a limited set of foreground spectra. 

The model for the binned CMB angular spectrum, C^ MB , is parameterized by cosmological parameters 
including the tensor-to-scalar ratio. Prior knowledge can be introduced in the modeling, for example by fixing 
the B-mode power spectrum of the lensed CMB to the known value. The noise spectra Njf are considered as 
additional components in each channel. Noise levels are either assumed known, in which case the noise spectra 
are fixed, or else are estimated simultanously with the foreground and CMB spectra. This model is fit to the 
measured angular spectra of multi-frequency sky maps. Since the spectra of each foreground component are 
expected to vary spatially, the fit is done in three spatial regions using apodized masks, to cut out the sharp 



4 The altICA algorithm is the result of a study supported by NASA LTSA Grant NNG04CG90G. 
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features of the Galactic plane. The model parameters, 9, are estimated by minimizing a mismatch statistic, 
(j), between the model and the data: 



E -f MC^r 1 ^) - logdet(C q (0)^C q ) 



(27) 



where C q (0) and C q are the binned spectra of the model and data in bin q, and w q is the effective number of 
modes in the bin. The errors on the CMB spectrum are estimated using the Fisher matrix. 

This method h as been used for CMB cleaning and power spectrum estimation for Planck simulations 
(|Leach et al.l 12008). It has also been used to estimate the temperatur e power spectrum of the CMB for 



WMAP first year data ( Patanchon et al. 2005 ) and for Archeops data ( Tristram et al. 2005h . To estimate 



the tensor-to-scalar ratio in the presence of foreground emission, one can assume a fixed shape of the CMB 
B-mode angular power spectrum, and then vary r in the parameter fit by varying the overall amplitude. This 
method produces estimates of the power spectra of each component. Maps of the CMB and of each foreground 
component can be estimated by Wiener filtering the observed multi-frequency maps, using the angular power 
spectrum of the best-fit model as the filter. 



3.4 Discussion 

These three classes of methods have all been applied to observed or simulated polarization data, and have 
different advantages and drawbacks. The template-cleaning method is simple and has been demonstrated 
to work effectively with current data. We do not yet have a good template for the polarized dust, or for 
low polarized components such as free-free, but multi-frequency information from CMBPol should be able to 
provide such templates. It is possible to propagate errors due to foreground uncertainty into the template- 
cleaned maps, but the framework for doing so is less well-defined than for some other methods. The parametric 
sampling methods are ideal for propagating errors due to foreground uncertainty in a fully Bayesian framework 
into both CMB maps and estimated power spectra. Additional observations can also be included as priors, 
such that all the information we have about the Galactic emission can be used. These methods could be 
biased if the physical model does not match the true sky, but various models can be tested to establish 
the most suitable one. They can also be numerically intensive. The blind component separation methods 
described here are quick to perform, and have the advantage of not assuming a particular physical source for 
the foregrounds. However neither do they take advantage of the fact that we do have useful knowledge of 
the emission processes, and can make unphysical assumptions about the foregrounds, for example including 
limited spatial variation of spectral indices. Cleaning currently observed polarization maps from WMAP has 
required methods operating in pixel space, so power spectrum methods may need further testing with real 
data. 

At the current time the community has tools ready to apply to polarization data with the sensitivity of 
a CMBPol satellite. Data from Planck and ground-based experiments will allow more extensive testing and 
comparisons between these methods, but we are ready now to tackle a CMBPol data set. In the event of 
a gravitational wave detection it is clear that having more than one method available will be important for 
cross-checking results. 



4 Numerical forecasts 

In this section we attempt to estimate whether a detection of a gravitational wave signal with r = 0.01 is 
possible for a realistic example mission. The specifications for this mission are given in Table [3] and have 
arisen from CMBPol studies focused on systematic effects and technological capabilities. We also consider 
how to optimize the allocation of detectors with frequency for such a mission. 
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Frequency 


FWHM 


AT (pol) 


GHz 


arcmin 


/iK / arcmin 


30 


26 


19.2 


45 


17 


8.3 


70 


11 


4.2 


100 


8 


3.2 


150 


5 


3.1 


220 


3.5 


4.8 


340 


2.3 


21.6 



Table 3: Specifications for a CMBPol 2m mission, for frequencies below 350 GHz. The noise limits 
AT are polarized limits for Q and U. 



4.1 Constraints on the CMB signal 

We consider two approaches to obtain forecasts on the B-mode CMB signal: a Fisher matrix method assuming 
that the foregrounds can be cleaned to a specified level, and a simulation approach, where we apply cleaning 
methods described in Section [3] to simulated maps, to estimate the CMB signal. For simplicity we will focus 
on the possible detection of a primordial signal with tensor-to-scalar ratio r — 0.01. Thi s level is picked 



as a ta rget that CMBPol should aim to reach, based on theoretical arguments discussed in iBaumann et al 



(|2008bl ). In this section all estimates are derived using the experimental specifications given in Table [3] 

Note that in addition to statistical and systematic uncertainties estimated in the approaches detailed 
below, the detection significance of a B-mode signal is also affected by cosmic variance. The error bars on the 
power spectrum, C#, with noise spectrum Ng, are given approximately by 



a(C e 



(21 + 1)/., 



sky 



-(C e + N e ). 



(28) 



As r tends to zero, the cosmic variance term also goes to zero. For non-zero r, cosmic variance will therefore 
affect the confidence interval. This means that a simple detection of B-modes can be confirmed with higher 
confidence than can a measurement of a specific small-amplitude signal such as r = 0.01. 



4.1.1 Fisher matrix approach 

The Fisher matrix approach described in Verde et al. ( 20061 ) is used for forecasting errors on measured cos- 
mological parameters. In this analysis the lensing contamination in the B-mode signal is accounted for by 
assuming that the lensing power spectrum can be accurately estimated, but not removed from the maps, 
such that it acts as an additional Gaussian noise term. To include foregrounds, the parameterization given 
in Section 12.51 is adopted for the total foreground signal, with the synchrotron given by Eqn. [6j The dust 
model is given by Eqn. [7] with Ad = 0.04 fiK 2 at I = 10, and a scale dependence m = —0.5, corresponding to 
an observed dust polarization fraction of ~4%. At scales I £ 30 this model matches the predicted spectrum 
from the Planck Sky Model. At larger scales the power spectrum of the simulated map flattens, which is not 
captured by this power law behavioi@. 



For a description of the method, see Appendix C in Baumann et al. ( 2008b ) and Verde et al. (j2006l ). 
The method assumes that at each frequency the angular power spectrum of the foregrounds can be cleaned 



5 This dust mo del corresponds to dust ' Model A' reported in lBaumann et al.l 12008b). An additional 'Model B' is 
also considered in IBaumann et alj 1 2008b ). with Ad = 0.002 fiK 2 at £ = 10, with scale dependence m = 0.6. The 



scale dependence of the dust is expected to lie within the —0.6 < m < 0.6 range considered in the two models. For 
simplicity we report only on Model A in this study, as it most closely matches simulated maps. 
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Method 



Average dust 
pol fraction (%) 



Description 



< 15 i < 150 



Fisher 
Fisher 





5 



No foregrounds 
10% residual 



0.0015 0.00046 
0.014 0.00052 



Parametric 
Parametric 
Parametric 



Fixed spectral indices 0.0015 
Power-law indices 0.0025 
Powerdaw indices 0.003 



0.0005 



Blind 



SMICA 



0.00055 



Table 4: The forecast la uncertainties on the tensor-to-scalar ratio for fiducial model r = 0.01, 
assuming other cosmological parameters are known, for a set of foreground assumptions. The 'Fisher' 
and 'Blind' tests include a residual noise signal from lensing. The 'Parametric' tests ignore lensing 
as a contaminant, a poor approximation at I > 15. For I < 15, the errors in the presence of 
foregrounds from the Fisher matrix test are higher than those obtained by other methods, as the 
assumed residual foregrounds are also higher at large scales. 



to a given residual level. The lowest and highest frequencies are discarded, assuming that they will act as 
foreground templates. Residual levels of both 30% and 10% of the total amplitude (10% and 1% in power) 
are considered, shown in Figure [T31 for a 10% residual at 90 GHz. These are chosen to represent a pessimistic 
and realistic case. Experience with previous experiments suggests that reducing the foreground amplitude by 
a factor of 10 (or more) should be possible, corresponding to the realistic case. The pessimistic case allows 
for poorer foreground cleaning, or a higher total foreground level than given by the model. 

Considering only the large-scale power, with I < 15, a model with r = 0.01 can be detected at almost 
7(7 in the absence of foregrounds, and assuming all other parameters are known. The error, a r — 0.0015, 
is given in Table [4] and is unmarginalized, so all other cosmological parameters are assumed known. With 
10% residual foreground amplitude, the primordial signal is undetectable according to this method, with 
o> = 0.014. However, this should be considered an overly pessimistic scenario, as improved cleaning is likely 
possible, and more importantly the large-scale foreground power is likely over-estimated, as the spectrum of 
the dust intensity flattens at large scales rather than following a pure power law. Including smaller scales, 
prospects are much improved. In the absence of foregrounds the r — 0.01 signal is detected at 20a with 
a r = 0.0004(H, and this is only decreased to 14<7 with pessimistic residual foregrounds. Wi th 10% residual 
foregro unds, the estimated degradation of the error on r is small. The analysis, reported in iBaumann et al.l 



(|2008bf ). also shows that information about the scale dependence of the tensors can be extracted by observing 
a range of scales, with forecasted errors on the tensor spectral index, n t , of a(nt) ~ 0.1. 

These forecasts should be interpreted with care. Some of the assumptions mad e in produci n g the se 



forecasts, for example, whether the lensing treatment is realistic (discussed further in ISmith et alJ ((2008)), 
and the assumed polarized dust levels at I — 100, need to be tested and may affect the resulting error bars. The 
estimated errors also assume that there is no effect of leakage of power from E to B modes. By using a large 
fraction of the sky, the errors on the measured polarization will vary spatially when foregroun d uncertainty 



i s incl uded, resulting in additional contamination of the B-mode signal from E-mode leakage. lAmarie et al 



(2005) approach this problem by masking out all regions of sky where the foregrounds exceed some fraction 



B This limit depends somewhat on assumptions made in the forecast. For example, in the absence of lensing noise, 
limits of ay = 0.0002 are obtained for f 3 k y = 0.87. 
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Figure 13: Comparison of CMB signals (E-mode - solid black; B-mode from inflation - solid blue,; 
total B-mode, including lensing - short dash), and an estimate of the residual foreground signal. 
Foreground spectra (for synchrotron plus dust) are given at 90 GHz for the cleanest 75% of the 
sky. The red line shows the total foreground signal, for simulated maps with mean dust polarization 
fraction 1.5%. Using a component separation technique on low resolution simulated maps gives an 
instrument noise plus foreground residual power spectrum (green), that is about 5% of the expected 
foreground at large scales. The green dash line is the power spectrum of foreground residuals 
assumed in the Fisher matrix analysis, where the mean dust polarization fraction is 5%. This likely 
over-estimates the residual foreground power at low t. 

of the CMB signal. This leads to larger sky cuts than the 20% cut used in this analysis, and larger errors 
than obtained here for the i < 150 case. We address this issue further in the next sub-section. 

4.1.2 Component separation with simulated maps 

An alternative approach is to apply foreground removal methods to simulated maps, in order to forecast the 
accuracy to which the CMB signal is recovered. In this initial study we test the parametric and blind com- 
ponent separation methods, as numerical codes already exist to do so. We do not consider template-cleaned 
forecasts here, leaving this approach for future investigation. 

Parametric methods: 

Simulations at each frequency for the CMBPol mission specified in Table [3] are formed using HEALPix with 
-^Vside = 16 and iV S idc = 128, using the following model: 

m(v) = A synch (u) + A dust (^) + f{v)K cxnh + f(u)n(v), (29) 

wherewhere f(v) converts thermodynamic to antenna temperature. The maps are formed using Eqn. [5] for 
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the synchrotron, and Eqn. [H]for the dust. Since they include the synchrotron and dust geometric suppression 
factors, the average dust polarization is 1% in the fiducial case (current PSM maps have average 5% polariza- 
tion fraction; we also consider maps with this increased dust level). Unlensed CMB maps are generated using 
the HEALPix synfast routine, for r = 0.01 with all other parame ters fixed to the best-fit ACDM parameters 
estimated from the WMAP five-year data Komatsu et al. ( 2008h . The N s id c — 16 (128) maps are smoothed 
with a 7 (1.5) degree beam. Gaussian noise is added uniformly over the sky, and 1// noise is ignored. In 
Figure [T3l we plot the angular power spectrum of the total foreground signal in these maps at 90 GHz. Part 
of the goal of performing the component separation is to assess the size of the errors on the estimated CMB 
maps, or power spectra, compared to the total signal. 

Marginalized CMB maps: We apply the parametric fitting methods described in Section [3.21 to the simu- 
lated maps, to estimate the foreground-marginalized errors on the observed CMB signal in map space. For a 
single random individual pixel outside the Galactic plane, the parameters of the model are sampled using the 
MCMC method described in Sec l3.2. Jf . and the marginalized error and bias on the CMB pixel amplitude are 
reported in Table El averaged over 5000 simulations. The errors are inflated about seven times compared to 
the foreground-free case. This level depends somewhat on the size of pixel in which the spectral indices are 
allowed to take unique values, in this case two degrees. We also estimate full marginalized CMB maps for 
this models The marginalized error map for the Q Stokes parameter is shown in Figure [14] In this case the 
spectral indices are fit in pixels of side 3.7 degrees, so errors are slightly reduced compared to the 2 degree 
pixels used in Table [5] The errors should be compared to the co-added noise level of 8 nK per pixel, in the 
absence of foregrounds. Increased errors, up to a maximum of 50 nK, are due to foregound marginalization. 
We then ask how big the marginalized errors on the CMB signal are compared to the total foreground signal. 
To do this we take the marginalized error maps and compute their angular power spectra, masking 25% of 
the sky using the WMAP P06 mask. This spectrum is shown in Figure Q21 and is comparable to the CMB 
spectrum for r = 0.01. The amplitude is about 5% of the total foreground signal at 90 GHz. This estimate 
may be over-optimistic if the true sky has more complicated spectral behavior and additional low polarization 
components that were not included in this test. However, in this test we allowed each pixel of side 3.7 degrees 
to have unique synchrotron and dust spectral index parameters. If we are only interested in a broad-band 
detection of the B-mode signal, larger pixels may be averaged together and could lead to lower residual error 
levels. 

Since there is about an order of magnitude uncertainty in the dust amplitude, we repeat the tests with 
varied foreground levels. For the case of an individual pixel, if we inflate the synchrotron amplitude three 
times, this inflates the errors by a further 50%, but a similar increase in the dust levels, where we are far less 
certain, has a smaller effect. If the simulation is modified to have a curvature in the spectral index, together 
with a 1% polarized free- free and 2% polarized spinning dust component with random polarization angles, the 
marginalized error is negligibly different from the fiducial case. In this case the model is sufficiently flexible 
that the 'synchrotron' component absorbs the additional polarized components, with an effective spectral 
index that is the weighted sum of the components. In future studies it will be useful to explore simulations 
that deviate further from the model, to test whether there is a significant bias in the recovered CMB power 
when the model fails to match the simulated data. 

To estimate the CMB power in these simulated maps, we use the Gib bs sampling method described in 
Section [3.2.2i using the Commander numerical code ( Eriksen et al. 200 7a ^1. Table [4] shows la errors on r, 
assuming other parameters are perfectly measured. When the foreground spectral indices are known perfectly, 
the tensor-to-scalar ratio is recovered with r = 0.0106 ±0.0017 for I < 15, and r = 0.0101 ±0.0005 for i < 150. 
These error bars are similar to those obtained in the Fisher matrix tests in the absence of foregrounds, where 

Using the FGFit code described in lEriksen et all ( 2006). 



Using the method and code described in lDunklev et al 



(20081 



We cross check results for low resolution maps by computing the exact likelihood of the marginalized maps, using 
the method in Section 13.2.11 These two methods should give the same results for any given set of simulated maps, 
although the implementation is different. 
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Case 


Error (nK) 


Case 


Error (nK) 


No foregrounds 


8 


Increased synchrotron 


79 


Fiducial model 


58 


Increased dust 


60 


Extended model 


58 


Increased dust + synchrotron 


81 



Table 5: Expected per-pixel sensitivities in polarization in 4x4 degree pixels in nK, including fore- 
ground degradation, for the 2m CMBPol example mission. The 'extended' model has a break in 
the synchrotron index, and low levels of polarized free-free and spinning dust emission. 'Increased' 
means a three-fold increase in amplitude. 



lcnsing noise is also included. These errors are too optimistic however, as in reality the synchrotron and dust 
spectral indices are uncertain and should be marginalized over. However, this provides a test of the method 
and shows that in this simple case the recovered signal is not significantly biased. It also indicates that regions 
of the sky wher e the foregrounds are high can still be used to estimate the CMB power: rather than masking 



out pixels as in lAmarie et alj (|2005f ). the errors are inflated in regions where foreground uncertainty is high, 
and propagated into the likelihood. 

A more realistic case includes the possible variation of synchrotron and dust spectral indices in each pixel. 
For pixels of side 4 degrees, and scales £ < 15, the error bar on r is increased to ay = 0.0025, indicating only 
a 4(T detection of the primordial signal. This is increased to oy = 0.003 when the dust component in the sim- 
ulated maps is made three times larger (corresponding to 4% polarization fraction). The limits may worsen 
if additional freedom in the foreground model is allowed (e.g. polarized free-free emission, a two-component 
dust power law, or a synchrotron spectral curvature), and if the dust is more then 5% polarized, although 
additional observations may provide priors or templates that would improve limits. Observations from the 
Planck satellite and ground based experiments will help to guide such choices for future preparatory studies. 
The application of this spectral index marginalization to the higher resolution maps, using £ < 150, was not 
calculated in this initial study, but is an obvious next step. 

Blind separation methods: 

The ILC method in spherical harmonic space described in Section [3. 3. II and the SMICA component separation 
method described in Section 13.3.21 have also been applied to simulated maps. 

For tests with SMICA, alternative maps to the fiducial simulation were used, as part of a Planck prepara- 
tory study. The results can be compared with others considered in this Section, as the maps are very similar 
to our 'Increased dust' model, with dust that is 4% polarized on average. The SMICA method assumes perfect 
knowledge of all cosmological parameters except r, the latter assumed to be measured only f rom B-modes 



For t he mission specifications given in Tabled limits of oy = 0.00055 are found for I < 150 (jBetoule et al 



20081 ) for a model with r = 0.01. This includes additional noise from the lensing signal, which when ignored 
leads to a limit of ay = 0.00046. Point sources are assumed to have a negligible impact on limits on r. These 
are similar to the limits found using the Fisher matrix method. When an additional 500 GHz channel is 
included, estimated limits improve to ay = 0.0003, indicating that this longer lever arm could provide useful 
additional information. However, it is possible that extrapolating from these high frequencies can introduce 
errors due to incorrect modeling assumptions. The numerical experiment with SMICA indicates that even 
with some complexity in the foreground modeling, estimation of B modes may be performed satisfactorily on 
a significant fraction of the sky, supporting large sky coverage for the CMBPol mission. 

Using the ILC method, rather optimistic levels of detection are obtained. Confidence levels on a model 
with r = 0.01 were not considered, but the 99% confidence level detection of r for fiducial r = is 0.0002 
for I < 8, and 0.0004 for 50 < £ < 130 in the absence of lensing0 This would imply foreground cleaning 



10 With lensing as a contaminant, this increases to r < 0.001 for £ < 8 and r < 0.008 for 50 < £ < 130. 
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Figure 14: Marginalized errors for a CMB Q Stokes map, using the CMBPol mission specified in 
Table [3] and the fiducial foreground model. In the absence of foregrounds the noise level would be 
8 nK in each pixel. The errors are increased where the signal is most uncertain due to foreground 
contamination. 

to better than ~1% in amplitude at large scales and better than 5% at intermediate scales, which is likely 
difficult to achieve in practice. 

4.2 Optimal use of focal plane 

Variables that drive experimental design include frequency range, bandwidth, and number of channels. To 
investigate the optimal configuration one requires either an analytic form for the errors, or a component 
separation code that can be run quickly. For this work, we employ the parametric method described in 
Section 13.2.11 to perform component separation for a single pixel, repeated many times with different CMB 
and noise realizations^ The averaged best-fit value for the CMB amplitude, and its marginalized error, are 
estimated. The error is then used as the statistic to compare the relative performance of different experimental 
specifications. In principle the number of channels, positioning in frequency space of the channels, and angular 
resolution can be varied. 

Here we perform a limited study for a fixed number of ten channels, chosen to allow for at least three 
parameters in each pixel for both synchrotron and thermal dust, and two for additional low-polarization 
components. The channels are spaced at logarithmic intervals in frequency, with the number of detectors at 
each channel varied such that the total signal-to-noise ratio ('total' is CMB plus foregrounds) is kept constant. 
The maximum and minimum frequencies are then varied, and the error on the CMB amplitude estimated. 
The foregrounds are assumed to have power law spectral indices for the parameterized model. The preferred 
minimum frequency is found to be at ~ 40 ±5 GHz, with a poorly constrained maximum at ~ 200 — 300 GHz. 
This spans a similar frequency range as proposed in the initial mission design, although has more channels. 
A number of different assumptions are investigated for the distribution of detectors at each frequency, using 
different scalings for sensitivities as a function of frequency They give similar results for the optimal frequency 
range. 

To use these tests to guide mission design, further investigation of the poorly constrained upper frequency 
limit would be useful. An extended parameterization for the foreground model could also be explored, to test 
how the optimal range is affected by the assumption of power law behavior. The number of channels can 
be explored with a similar tool, but was not considered in this initial study. It is clear that more channels 
will help to remove the foregrounds effectively and break degeneracies between spectra of different emission 
mechanisms. Sparsely placed frequency channels lead to uncertainty in the curvature of the foreground 
frequency spectra. However, there is always a trade-off between increased sensitivity at each frequency, and 
increased number of channels. 

11 FGFIT_pix was designed exactly for this purpose, with input files for multiple experimental design and sky models. 
The experiment file contains the frequency and noise level for each channel. 
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If we ignore the effects of lensing there is little motivation to go to higher resolution than 1 degree. 
However, optimization for lensing should be considered in order to extract the gravitational wave signal at 
smaller s cales where the for egrounds are lower than at the lowest multipoles. See the companion lensing 
report bv lSmith et al.l (2008) for more details. 



4.3 Future tests 

The forecasts presented here represent an initial investigation into the CMBPol mission capabilities. Work 
remains to be done to make concrete forecasts, but we will be in a better position as more data is obtained 
from the Planck satellite and ground-based experiments. Further preparatory work for a CMBPol mission 
could include: 

(i) Bias and modeling: So far we have examined cases where the simulated maps are correctly described 
by the models. In a realistic scenario, the model will not perfectly capture the observed two-dimensional 
emission. It will be useful to determine how many parameters we need to fit the true sky 'well enough', and 
whether an incorrect model leads to significant bias in the estimate of the CMB signal. 

(ii) Convergence of estimates: Different methods exist for forecasting the constraining abilities of the 
CMBPol mission. Consensus of the effect of foregrounds on the estimated B-mode signal is important in the 
planning of the mission, and will likely improve with better observations. However, not all methods will clean 
the maps equally well, so we would not necessarily expect to achieve the same errors on the CMB signal from 
all component separation methods. We plan to continue with the program we have begun, to compare the 
estimated residual noise-plus-foreground levels for various methods over the I < 150 range, accounting for the 
lensing signal as a source of noise in simulations. These would include the Fisher methods, parametric and 
blind separation methods described here, as well as template cleaning applied to forecasting. 

(iii) Optimizing the experiment: The community has tools that can be applied to experimental opti- 
mization. In particular the parametric fitting employed in the 'FGfit' code allows for easy comparison of the 
CMB errors in a pixel, marginalized over foregrounds, for different experimental set-ups. Other methods, 
including template cleaning and ILC component separation, can be used for alternative investigations. Fur- 
ther exploration with these tools will be useful, including tests for the optimal number of channels, as well as 
variation of the complexity of the model. One can also test the effect of varying the experimental sensitivity, 
to determine the noise level at which the observations become foreground dominated. 



5 Summary 

Polarized Galactic emission at microwave frequencies is dominated by two components, synchrotron emission 
and thermal dust emission. The synchrotron emission has been well measured on large angular scales at 
23 GHz by WMAP. Our knowledge of polarized dust emission is relatively poor, particularly in the low 
surface brightness regions out of the Galactic plane targeted for CMB observations. It will soon be better 
characterized by the Planck satellite, due for launch 2009, which will measure polarized emission with a spatial 
resolution < 30 arcmin at 30 - 350 GHz. With the assumption of a polarized dust fraction of 1-10%, the total 
foreground emission is expected to be roughly a factor of eight brighter, on angular scales of two degrees over 
75% of the sky, than a gravitational wave signal with tensor-to-scalar ratio r = 0.01. There may also be small 
polarization contributions from free-free emission, spinning dust emission, or other 'anomalous' components 
that have not yet been detected. 

Past experience working on analogous problems suggests that the challenge of dealing with significant fore- 
ground components is not insurmountable, and has been addressed successfully with previous experiments 
such as FIRAS. In this case the CIB foregrounds were reduced by a factor of ten to extract the signal, com- 
parable to the cleaning required of polarized Galactic foregrounds for the CMB. Although past performance 
is no guarantee of future success, we have consistently been able to clean foregrounds to within a factor of 
a few times the uncertainties of the raw measurements. As a community we have a set of tools ready for 
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performing component separation and estimating the CMB signal from observed sky maps. These could be 
applied today to CMBPol data, and will be significantly improved and tested with data from Planck and 
from ground and balloon-based experiments during the next five years. There is some variation in predicted 
performance for a specified CMBPol mission, so the refinement process with upcoming experiments will be 
useful for understanding why some methods perform better than others, and whether this depends on details 
in our simulations. 

Despite the good prospects for cleaning foregrounds, such prospects rely sensitively on experimental design. 
Observations over a wide range of angular scales should be planned in order to mitigate anticipated residual 
Galactic contamination on the largest angular scales (I < 15). Multiple channels are needed to characterize 
and remove the foregrounds. Although there are only two dominant polarized components, they are unlikely 
to have spectral indices that can be modeled accurately as power laws. This means that at least three channels 
per component should be included. Since there may also be additional low-polarization components, at least 
ten channels spaced logarithmically in the range ~ 40 < v < 200 — 300 GHz would be advantageous, although 
further studies should be undertaken to explore this in more detail. Extrapolating from very low and high 
frequencies can introduce modeling errors, so frequencies below ~ 300 GHz are more easily used for foreground 
subtraction. 

Current methods and simulations indicate that we can clean Galactic foregrounds from maps of the po- 
larized sky to at least the 5 — 10% level. There are therefore good prospects that, for a realistic CMBPol 
mission design, a gravitational wave signal with r = 0.01 could be detected at more than 10er. Successfully 
disentangling this signal from the emission of our own Galaxy will provide rich rewards, as a detection of 
primordial B-modes would open up a new window on the earliest moments in the universe. 
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